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Abstract 

In this article, we study the three-pion B c -decays £?+ — > J/ij) n + n~n + , r) c n + n~ n + with 
dominance of the intermediate axial-vector meson ai(1260) and vector meson p(770) in the 
tv + -k~7v + and tv + -k~ invariant mass distributions respectively, and make predictions for the 
. branching fractions and differential decay widths. The ratio between the decays — > 

J/ip 7r + 7r~ 7r + and — > J/ip7v + is compatible with the experimental data within uncer- 
tainties, other predictions can be confronted with the experimental data in the future at the 
LHCb. 
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^ ! 1 Introduction 

In 1998, the CDF collaboration observed the pseudoscalar Bf mesons through the semileptonic 
decays Bf — > J/tpl^X and Bf — > J/ipi^Dt in pp collisions at the energy y/s = 1.8 TeV at the 
Fermilab Tevatron pQ. The bottom-charm quarkonium states B c are of special interesting, the 
ground states and the excited states lying below the BD, BD*, B*D, B*D* thresholds cannot 
annihilate into gluons, and therefore are more stable than the corresponding heavy quarkonium 
states consist of the same flavor, and would have widths less than a hundred KeV [2j. The 
semileptonic decays Bf — > J/t/j^Dg, B/ — > J/ipe + D e were used to measure the B c lifetime [21 S], 
which is about a third as long as that of the B and B s mesons as both the b and c quarks 
ly-j | decay weakly. In 2007, the CDF collaboration observed the B^r mesons with an significance 

■ exceeds 8cr through the decays Bf — > J/ipir in pp collisions at the energy y/s — 1.96 TeV, and 

obtained the value m_B c = (6275.6 ±2.9 ±2.5) MeV [5]. In 2008, the DO collaboration reconstructed 
{SJ ■ the decay modes Bf^ — > J/ipir^ and observed the B^r mesons with an significance larger than 

5cr, and obtained the value = (6300 ± 14 ± 5) MeV [6]. Now the average value is = 
(6.277 ± 0.006) GeV from the Particle Data Group [7 . Recently, the LHCb collaboration observed 
the decay B/ — > J/-07r + 7r~7r + for the first time using 0.8fb _1 of the pp collisions at y/s = 7 TeV, 
the measured ratio of branching fractions is 



Br(B+ -> J lbTT + 7T-1T + ) 

— v c t , 7—^ = 2.41 ±0.30 ±0.33, (1 

Br (Be -> J/ipn+) V ' 

where the first uncertainty is statistical and the second is systematic [8]. 

The hadronic decays B/ — > J/tp ir + Tr~ir + , J/ipir + take place through the weak 6-quark decays 
b — > cW* — > cud, which are analogous with the r-lepton decays t — > v T W* — > v T ud. We can 
use the existing experimental data on the r-lepton decays to obtain reliable prediction for the 
B/ — > J/ipir + Tr~ir + branching fraction with the method of spectral functions, 

d$(W* -> 7T + 7r-7r+)e AI e: = (q^ - q 2 g^)p T (q 2 ) + q^updq 2 ) , (2) 

where the d§(W* — > tt + tt^tt + ) is the Lorentz-invariant three-body phase factor, the e M is the 
effective polarization vector of the virtual H^-boson, the spectral functions pr(q 2 ) and Pl(q 2 ) are 
universal and can be determined by theoretical and experimental analysis of the r — > i/ T Tr + ir~ tt + 
decays [HI EH [HI H2] ■ The spectral function ph{q 2 ) is negligible according to conservation of vector 
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current and partial conservation of axial- vector current. The explicit expression of the spectral func- 
tion prig 2 ) can be fitted to the experimental data or calculated by some phenomenological models, 
the spectral function pr(q 2 ) is always saturated by exchange of the intermediate axial- vector meson 
ai(1260) with an special ansatz for the vertexes aipn and p7T7T based on some phenomenological 
quark models jT3l [14l USl UB] . In this article, we intend to study the decays — > J / ipir + ir~ ir + , 
77 c 7r + 7r~7r + , J/ipTT + , r) c Tt + with the phenomenological Lagrangians, and calculate the Feynman di- 
agrams directly. The decays Bj — > i] c ir + n~n + , t] c tt + have not been observed yet, but they are 
expected to be observed in the future at the Large Hadron Collider (LHC). The B c decays will be 
studied and their branching fractions will be determined in the LHCb experiments. 

In Ref.[17j. Lichard and Juran perform detailed analysis of the vertex aipir with the following 
Lagrangian, 

C = ^f- {cos8 [(d^p " - d"p 0tl )a~d u n + ] - sin6 [(d"p 0v - d v p^d^^] +■■■}, (3) 

considering the decays ai(1260) — > pit, where the momenta of the p and tt mesons in ccnter-of- 
mass frame of the initial ai(1260) meson are about 0.37 GeV. Such a Lagrangian maybe lead to 
amplified amplitude for the decay B+ — > J/ij) n + n~n + as the sub-amplitude a^~(1260) — > pn + can 
be accounted as g af mlp • k w or g af mlp • Qat where the momenta q a , l p and k v are large, we have 
to introduce form-factors to parameterize the off-shell effects. In the decays r — > ai(1260) v T and 
£?+ — » J/ip af (1260), the momenta of the ai(1260) meson in center-of-mass frame of the initial 
particles are about 0.45 GeV and 2.16 GeV, respectively. On the other hand, we know that the 
decays ai(1260) — ¥ pir are S-wave dominated [7], we prefer the simple Lagrangian, 

£a P 7r = gap7ra1Pp.Tr, (4) 

in this article. Furthermore, we use the Lagrangian, 

CpTTTT = -igpTriTp at2 Tr~(£) p - ti)^ + , (5) 

to study the vertex pirir [18] . 

The article is arranged as follows: we derive the decay widths of the processes B+ —> J/tp tt + it~it + , 
rj c iT + Tr~7T + , J/'07r + , i] c 7r + in Sect.2; in Sect. 3, we present the numerical results and discussions; 
and Sect. 4 is reserved for our conclusions. 

2 Decay widths of the processes — > J/i/jtt + 7t ~tt + , r] c 7T + 7T^7T + 

The hadronic decays B+ — > J/iJjtt + , t] c tt + , J/ijj tt + tt~tt + and rj c n + ir~n + can be described by the 
effective Hamiltonian, 

n cS = ^VcbV^dip^il-^bdni-^u + h.c, (6) 

where the V c b, V u d are the CKM matrix elements, the Gf is the Fermi constant, and the C\(p) 
is the Wilson coefficient defined at an special energy scale, Ci(mb) = 1.14 [19 . In the following, 
we write down the definitions for the weak form-factors Fi(q 2 ), Fo(q 2 ), Ai(q 2 ), ^(q 2 ), A?,(q 2 ), 
A (q 2 ) and V(q 2 ) for the current J p = cj^(l — -f 5 )b sandwiched between the B c and the rj c , J/tp 
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( Vc (p)\J^O)\B c (P)) 

m(p)\j^o)\B c (p)) 



In 

,J J/ip) q 2 



(q 2 ) + (rnk-m 2 j^) ^F (q 2 ) , 



(P+ P )^F + (q 2 )+q^(q 2 ), 

% |-e£(m Be + mj/^A^q 2 ) + e* ■ P(P + p\ 



(7) 



A 2 (q 2 



2mj^e* ■ P% [A 3 (q 2 ) - A (q 2 )] - ie^e* v q a (P + p) 



v(q 2 



e;(m Bc + mj/^A^q 2 ) + e* ■ P(P + p)„ 
A_(q 2 ) 



m Bc + mj/Tp 



ie^e* u q a (P + p) 



A+(q 2 ) 
m Ba + mj/,p 



(8) 



where 



F+(q 2 ) = 

F-(i 2 ) = 

A 3 (q 2 ) = 

A + {q 2 ) = 

A-(q 2 ) = 



Fl(q 2 ), 

[Foiq^-F^q 2 )] {m% c 



m Bc + mj/^ 



Al{q 2 ) 



m Bc - mj/j, 



2mj/^ ' 2mj/,p 



A 2 (q 2 



A2{q 2 ) , 

2m Jf ,( mB;+mj/ ,) 



(9) 



and Vb(0) = V^O), and the e M is the polarization vector of the vector meson J /if). 

There have been several approaches to calculate those weak form-factors, such as the QCD sum 
rules PI1EUE2], the light-cone QCD sum rules [23], the relativistic quark models [24 ] 125 ] 126 ] [27] . 
the light- front quark models [SHIES], the non- relativistic quark models [30] [31], perturbative QCD 
[32] . etc. In this article, we take the typical values from the QCD sum rules [19], the relativistic 
quark models [Ml [213 an d the light- front quark models [55], and refer them as QCDSR, RQM1, 
RQM2 and LFQM, respectively. In calculations, we use the following definition for the decay 
constant f a of the axial- vector meson ai(1260), 



|u(0)7^75<0)|a 1 (1260)) = f a m a e 



n ■ 



(10) 



where the is the polarization vector. 

We take into account the effective Hamiltonian H c s, the weak form- factors and the Lagrangians 



£-ap-n, C-p-K-n to obtain the amplitudes Tj/^^+j,.- ,,.+ , X , 1)c7r + 7r - 7r + of the processes — ¥ J/ipir + ir n + , 



r/ c 7T T 7r 7T T , 



T 



Gp 



V2 

fj = T] + T 2 , 



V c6 T/ [1 * d Ci(mb)/ a m a r j , 



(11) 



3 



where j — J/tjj tt + tt 7r + , ?7 c 7r + 7r tt + , 



-W) , ,* Pfl -Mg 2 ) 



»,x,..< Vi /' - Pf V[ f ) ) ~T=z~TT7 (-<T + ) <</ , 



+ ) iff^Ci - r) a , (12) 



/ 2 -m 2 p +iVPT p (l 2 ) V I 
rj/^-^ = T^ 7r+7r _ 7r+ ( 7 r+(r)0 7r+(fc)) 1 (13) 



q 2 - m 2 + iy/q 2 T a (q 2 ) \ ' q 



Tl^ + = {(P+p),F + (q 2 )+q,F_(q 2 )}- | - ( -g^ + ^ ) ^ 



S^a + ) ig^{t- r) Q , (14) 



^+.-.+ = ^-^(t+W^t+W), (15) 



then obtain the differential decay widths 

IT- 1 2 da 2 dl 2 

dT j = ±^^-?Ld$(P^ q>P )d$(q^l,k)d$(l^r,t), (17) 



where the d$(P — > q,p), d$>(q — > i, k), d$(Z — >• r, t) are the two-body phase factors defined 
analogously, for example, 

d<£(P^q,p) = (2*)H\P-q-p) % . (18) 

(2tt)' ;, 2po (27r) d 2(?o 

The decay widths of the processes 2?+ — > J/tpir + , r/ c 7r + can be calculated straightforward using 
the effective Hamiltonian and the weak form- factors, the explicit expressions are neglected for 
simplicity. 



3 Numerical results and discussions 

The input parameters are taken as G F = 1.166364 x l(T 5 GeV~ 2 , V ud = 0.97425, V cb = 40.6 x 
10" 3 , m T = 139.57 MeV, m p = 775.49 MeV, T p = 146.2 MeV, m Ba = 6.277 GeV, r Ba = 0.45 x 
10~ 12 s from the Particle Data Group [7], m a = 1255 MeV, T a = 367MeV from the COMPASS 
collaboration [33], = 6.05 from the decay p —> irn 18 J, g ap ix = 3.37 from the decay ai(1260) — > 
pit 0, f a = 0.24 GeV from the QCD sum rules [34]. 

We obtain the branching fractions of the decays B£ — > J/t/jir + , ?/ c 7r + , J/tp tt + tt^tt + and 
?; c 7r + 7r _ 7r + with the typical values of the weak form-factors from the QCD sum rules [19], the 
relativistic quark models [231 HS]i an d the light-front quark models [28 . The form- factors in 
Refs. [HI [Ml 113 are fitted to an single pole form, 

/(9 2 ) = -1 f ~§ 2 , (19) 
1 - q 2 /mj it 
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0.63 
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0.66 
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[4 51 


[4 51 


[4 51 


[4 51 


RQM1 [24] 


0.56 


0.54 


-0.95 


0.83 


0.61 


-0.32 




[5.45] 


[4.76] 


[4.68] 


[4.72] 


[4.84] 


[4.80] 




Ax(0) 


4,(0) 


A Q (0) 


V(0) 


Fi(0) 


F o (0) 
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[m /it ] 


[m/it] 


[m/it] 


[m f u] 


[m f i t ] 


RQM2 [25] 


0.50 


0.73 


0.40 


0.49 


0.47 


0.47 




[4.84] 


[4.72] 


[4.04] 


[3.99] 


[4.41] 


[4.72] 




MO) 


A 2 (0) 


A Q (0) 


no) 


Fi(0) 


F o (0) 




[ci/c 2 ] 


[ci/c 2 ] 


[ci/c 2 ] 


[ci/c 2 ] 


[ci/c 2 ] 


[ci/c 2 ] 


LFQM [IS] 


0.50 


0.44 


0.53 


0.74 


0.61 


0.61 




[1.73/0.33] 


[2.22/0.45] 


[2.39/0.50] 


[2.46/0.56] 


[1.99/0.44] 


[1.18/0.17] 



Table 1: The parameters for the weak form- factors, where the unit of the irifn is GeV. 



while the form- factors in Ref . |28j are fitted to an exponential form, 



f(q 2 ) = f(0)exp( Cl q 2 +C 2 q 4 ) 



(20) 



where the f(q 2 ) denote the weak form-factors, the m/it, ci, c 2 are fitted parameters. The numerical 
values are presented explicitly in Table 1. 

The numerical values of the branching fractions are shown in Table 2, from the table, we can see 
that the branching fractions vary in a large range according to the values of the weak form-factors 
from different theoretical approaches, it is difficult to determine which one is superior to others. 
In Table 2, we also present the predictions from the Berezhnoy-Likhoded-Luchinsky (BLL) model 
for comparison, where the spectral function 



Pr{q 2 



x 10" 



1 



9mz 



1 + 190<7 2 



(21) 



determined from the r 



V r Tr + TT 



[{q 2 - 1.06) 2 + 0.48]^ 

r + decays is used (9[ [TUl E] ■ The present values are slightly 
different from that of Ref. [5], as we have taken slightly different input parameters. 

The ratios among the branching fractions are shown explicitly in Table 3, from the table, we 
can see that 



Br{B+ -> J/V>7r+7i--7r+) 
Br{Bt ->• J/V>7r+) 



2.17, 1.87, 2.28, 1. 



(22) 



are all compatible with the experimental data 2.41 ± 0.30 ± 0.33 within uncertainties jS], while 
the ratios based on the weak form- factors from the QCD sum rules in Ref.piJ] and the relativistic 
quark models in Ref. [2 5) are better. All those predictions can be confronted with the experimental 
data in the future at the LHCb collaboration. 

In Fig.l, we plot the differential decay widths of the B c mesons dT(B+ — > J/ipir + ir~ir + )/dq 2 , 
dT{B+ -> r lc n+Tr-Tr+)/dq 2 , dT{B+ -> J/<ip tt+tt- tt+) / dl 2 and dT(B+ -> Vc n+ir- tt+) / dl 2 with 
variations of the squared momenta q 2 (M 2 (7r + 7r _ 7r + )) and I 2 (M 2 (tt + n~ )) . By measuring the 
momenta dependence of the differential decay widths, we can test dominance of the intermediate 
axial- vector meson ai(1260) and vector meson p(770) in the invariant 7r + 7r~7r + and 7r + 7r~ mass 
distributions, respectively. It is difficult to distinguish the two tt + mesons in the final states, so we 
should not be serious for taking the squared virtual momentum I 2 as the invariant mass distribution 
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QCDSR .19 


RQM1 .24 


RQM2 25, 


LFQM [28 


B+ -> J/^7T+ 


1.222 


1.106 


0.544 


0.956 


B+ -> rjc 7T+ 


1.592 


1.360 


0.879 


1.360 




2.655 [3.475] 


2.072 [2.758] 


1.242 [1.607] 


1.797 [2.390] 


£>+ — > ry c 7r + 7r _ 7r + 


1.854 


1.544 


0.947 


1.578 



Table 2: The branching fractions of the B c decays, where the unit is 10 -3 , and the references 
denote the hadronic form-factors from that articles are used. The values in the bracket denote 
the predictions from the BLL model, where the spectral function pt{<] 2 ) determined from the 
t — > i> T TT + Tr~TT + decays is used. 





QCDSR 19 


RQM1 [24 


RQM2 [25, 


LFQM 28] 


B7'(B+^.//i/>7r + 7T~7r + ) 
Br(B+^J/ipTT+) 


2.17 


1.87 


2.28 


1.88 


Br(B+->r) c 7r+7r~7r+) 
Br(B+^J/ipTT+) 


1.52 


1.40 


1.74 


1.65 


Br(B+^jj c 7T + ) 
Br(B^^J/ifm+) 


1.30 


1.23 


1.62 


1.42 


Br(B+->r) c 7r+7r~7r+) 


1.16 


1.14 


1.08 


1.16 



Table 3: The ratios among the branching fractions of the B c decays, where the references denote 
the hadronic form-factors from that articles are used. 

M 2 (ir + Tr~). However, from Fig. 2 we can see that the approximation I 2 — M 2 (ir + ir~) works rather 
well. 

In Fig. 2, we plot the 7r + 7r~7r + and 7r + 7r~ invariant mass distributions in the decays B c — > 
J/?/!7r + 7r~7r + compared with the experimental data and the predictions of the BLL model. With 
suitable normalization, the weak form- factors from the QCD sum rules |19) . the relativistic quark 
models [5U [25] , and the light-front quark models 28 lead to almost the same line-shapes for the 
invariant mass distributions, although they correspond to quite different decay widths. From the 
figure, we can see that the present work and the BLL model both describe the experimental data 
on the invariant mass distributions M(n + ir~Tr + ) well, while the present prediction is better. 

4 Conclusion 

In this article, we study the three-pion i3 c -decays B£ — > J/tp tt + tt~tt + , 7/ c 7r + 7r~7r + assuming 
dominance of the intermediate axial-vector meson ai(1260) and vector meson p(770) in the invari- 
ant 7r + 7r~7r + and 7r + 7r~ mass distributions respectively, and make predictions for the branching 
fractions and differential decay widths. The ratios between the decays B^ — > J / 'ip 7r + it~ tt + and 
B+ — > J/ipir + based on the form-factors from different theoretical approaches are compatible with 
the experimental data within uncertainties. The decays B£ — > rj c tt + tt~tt + , r) c Tr + have not been 
observed yet, the predictions can be confronted with the experimental data in the future at the 
LHCb. 
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Figure 1: The differential decay widths of the B c meson, where the A, B, C and D denote 
the differential decay widths dT{B+ -> J/4>ir+TT-ir+)/dq 2 , dT{B+ -> t) c tx+ v+) j dq 2 , dT{B+ -> 
J /ipir + ii~ n + ) / dl 2 and dT(B£ — > rj c ir + ir~ir + )/dl 2 , respectively. 




Figure 2: The invariant mass distributions of the 7r + 7r tt + and 7r + 7r compared with the experi- 
mental data, where we have taken the approximation I 2 = Af 2 (7r + 7r~). 
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